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1.0 SUMMARY 


Program L225 (TEV126) is structured as three overlays, one main and two primary. 
Input into the program is made via cards and magnetic files (tapes or disks). Output 
from the program consists of printed results and magnetic files containing vectors 
suitable for plotting. 

Although L225 serves as a module of the DYLOFLEX system, it can be operated as a 
standalone program. Subroutines used by L225 include routines embedded in the 
program code, routines obtained from the standard FORTRAN library, and routines 
obtained from the DYLOFLEX library. 

2.0 INTRODUCTION 

The computer program L225 (TEV126) can be used as either a standalone program or as 
a module of a program system called DYLOFLEX (see fig. 1) which was developed for 
NASA under contract NASl-13918 (ref. 1). Because of the DYLOFLEX contract 
requirements developed in reference 2, a program was needed to calculate the time 
histories of different airplane responses to discrete excitations using frequency response 
data calculated in a power spectral density solution. An existing program^ was revised 
according to DYLOFLEX specifications'^ to make it compatible with other modules of 
the DYLOFLEX system. 

The objective of this volume is to aid those persons who will maintain and/or modify the 
program in the future. To meet this objective, the following items are defined in detail: 

• Program design and structure 

• Overlay purpose and description 

• Input, output, and internal data base descriptions 

• Test cases used in program checkout 


^Smith, B. E.; and Johnson, M. R.: Discrete Gust Convolution Program (TEV126). Boeing 
document D6-29668TN, January 1969. 

^Every, J. S.: Frequency Response Function Sorting Program (TEV126, Linkl , DIGl). Boeing 
document D6-29669TN, vol. II, April 1971. 

^Clemmons, R. E.: Programming Specifications for Modules of the Dynamic Loads System to 
Interface with FLEXSTAB. NASA contract NASl-13918, BCS-G0701, September 1975. 
(Internal Document) 




Figure 1. - DYLOFLEX Flow Chart 









3.0 PROGRAM DESIGN AND STRUCTURE 


The program is structured as a system of three overlays, one main overlay and two 
primary (see fig. 2). 

Main overlay (L225,0,0) L225vc 

Primary overlay (L225,l,0) SORT 

Primary overlay (L225,2,0) ANAL 

The main overlay L225 reads data cards that direct the execution of the primary 
overlays. It also aids communication between the primary overlays via labelled common 
blocks. 

The 1,0 primary overlay SORT reads the frequency response functions from cards or 
tape, sorts those functions the user wishes to use in the time history convolution, and 
saves the desired response functions on a temporary scratch file. 

The 2,0 primary overlay ANAL reads the sorted frequency response functions from the 
scratch file generated by SORT. The retained functions are then convoluted with a 
forcing function chosen by the user. Output can be in the form of printed time histories, 
printed maximum and minimums, or tapes containing data vectors suitable for plotting. 



Sort frequency 
response functions 
and retain those 
requested for 
analysis in ANAL 


Convolute the 
frequency response 
functions and 
prepare the 
time histories 


Figure 2.— Overlay Structure of L225 (TEV126) 
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The input and output for each of the overlays is shown in figure 3. 

Card Sets 

Card Sets 1*0, 1.1, 1.2, 2.0, Card Sets 



(PLTRES) 


Figure 3.— Input/Output of L225 Overlays 


Although L225 (TEV126) serves as a module of the DYLOFLEX system, it can be 
operated as a standalone program. When the program is run by itself, it becomes the 
user's responsibility to generate input data in the format required by L225 (TEV126). 
(See volume I of this document for the required data and formats.) 
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This program requires subroutines that are not part of the L225 code. Some routines are 
automatically obtained from the standard FORTRAN library when the program is 
loaded. Others, however, are stored in the DYLOFLEX alternate subroutine library 
which must be declared at the time of loading. Subsequent sections describe each 
overlay separately and contain tables displaying the routines called and the library in 
which they are located. Each subroutine has only one entry point. 

This volume describes the program in a macro sense. A more detailed discussion 
appears in the comments contained in the program source code. Each routine contains a 
preface describing the routine’s purpose, author, analytical steps, modification history, 
input data, output data, glossary of variables, and a list of other routines called. 
Embedded within the executable code are comments labeling each section and 
explaining logical branches. 


3.1 OVERLAY (L225,0,0) - L225vc 

The main overlay of L225 (TEV126) is L225vc, where v is a letter indicating the 
program version, and c is an integer number indicating the correction number which 
applies to the v version. 


Purpose of L225vc 

L225vc performs certain bookkeeping tasks, directs the execution of the primary 
overlays, and aids communication between primary overlays via labeled common blocks. 

Analytical Steps of L225vc 

L225vc performs its task in the following six steps: 

1. The subroutine PRGBEG is called to place the program header on the printed 
output. 

2. A data card is read. It must begin with $TIME to ensure that the card input file is 
correctly positioned. If it does not contain $TIME, execution is terminated. 

3 A program directive card is read, printed, interpreted, and acted on according to 
the following logic; 

If the keyword is $TITLE, jump to step 3 again. 

If the keyword is $CHECkout, set the checkout switch and jump to step 3 again. 

If the keyword is $SORT, jump to step 4. 

If the keyword is $ANAL, jump to step 5. 

If the keyword is $QUIT jump to step 6. 
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4. Overlay (L225,l,0) is called. When it is finished, program control returns to step 3. 

5. Overlay (L225,2,0) is called. When it is finished, program control returns to step 3. 

6. Subroutine PRGEND is called. This subroutine places the program termination 
message on the printed output. 

If fatal errors are discovered, the program prints a diagnostic and jumps to step 6. 

The macro flow chart of this overlay is shown in figure 4. The subroutines called are 
displayed in table 1. 


Table 1,— Routines Called by L225vc 


OVERLAY (L225,0,0) 
PROGRAM L225VC 


ANAL - OVERLAY (L225,2,0) 
FETAD + 

FETDEL + 

IRDCRD + 

PRGBEG + 

PRGEND + 

RETURNF + 

SORT - OVERLAY (L225,1,0) 


+ indicates a routine on the DYLOFLEX alternate subroutine 
library DYLIB. 

All others are local to L225 (TEV126) . 
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Figure 4.— Macro Flow Chart of Overlay (L225,0,0) L225vc 













I/O Devices of L225vc 


L225vc reads program directive cards and writes them on the printed output file. All 
other I/O accomplished by L225 (TEV126) is done in the primary overlays. 

3.2 OVERLAY (L225,l,0) - SORT 

Purpose of SORT 

SORT is the first primary overlay in program L225 (TEV126). SORT reads frequency 
response function data from either cards or a magnetic file (tape or disk), selects 
responses requested via card input, and writes a file containing the response functions 
to be analyzed in ANAL, the (L225,2,0) overlay. 

Analytical Steps of SORT 

SORT performs its task in the following steps: 

1. Subroutine RDSZSC is called to read the sorting instruction cards (card sets 3.0 
through 6.0). 

2. Subroutine CARDIN or TAPEIN is then called to read the first set of responses 
from cards or tape, respectively. 

3. Subroutine WRITAP is called to write the selected reponses on the scratch file 
SRTFRP. 

If the input response functions are on cards, control is shifted to step 6. 

4. Subroutine TAPEIN is called again to read the second set of responses from tape. 

5. WRITAP is called to write the second set of selected responses on SRTFRP. 

6. Control is returned to the calling program, the main overlay L225vc. 

When input is made via the standard input file NTAPE (see file map shown in vol. I), 
the input response functions (generalized coordinates and loads) are processed in two 
separate passes. Steps 2 and 3, the generalized coordinates (Q’s) are read, sorted, and 
written on the scratch file SRTFRP Steps 4 and 5, the load frequency response 
functions are read, sorted, and written on SRTFRP. This procedure is used because it 
saves core storage and allows the user to specify Q's and loads separately. 

The responses on NTAPE are defined for one frequency at a time; i.e., the tape consists 
of a set of arrays containing all the responses (both Q’s and loads) at frequency one, 
followed by a similar set of arrays at frequency two, and so on until the final set of 
arrays at frequency NFREQ. With the first call of TAPEIN, the real and imaginary 
parts of the generalized coordinate arrays are stored in the rows of the RESR and RESI 
matrices in SORT. These matrices are of size NFREQ by NRMAX, where NRMAX is the 
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larger of the two numbers, the number of generalized coordinates, or the number of 
loads. It then becomes the task of WRITAP to select the columns of RESR and RESI to 
be written onto SRTFRP for analysis in ANAL. The second call to TAPEIN repeats this 
process for the load frequency response functions. 

If the reponse functions are input via cards, the subroutine CARDIN is only called once. 
Response functions are read from cards and stored in successive arrays of size 1 by 
NFREQ. In this case, NRMAX is the sum of the number of responses analyzed (see card 
input for SORT in vol. I). No sorting is required with card input. The responses are 
written directly on the scratch file SRTFRP. 

Figure 5 shows the macro flow chart of the SORT overlay. The subroutines called are 
displayed in table 2. 


Table 2.— Routines Called by SORT 


OVERLAY (L225, 1,0) 

PROGRAM SORT 

CARDIN 

FNDKEY + 

LOCF * 

RDSZSC 

REQFL + 

TAPEIN 
WRITAP 

* indicates a routine on the FORTRAN subroutine library. 

+ indicates a routine on the DYLOFLEX alternate subroutine 
library DYLIB. 

All others are local to L225 (TEV126) . 


|lRDCRD + 

( IRDCRD + 
NAMFIL + 


FETAD + 

« FETDEL + 
READTP + 

WRTETP + 
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START 



Routine to read from 
NTAPE the first set of 
responses to be sorted 



NTAPE 


< 3 - 


Write the first set of 
responses on SRTFRP to 
be read and analyzed 
by ANAL 


Routine to read cards designating 
the problem size, data source 
and responses to be saved 



Routine to request field length 
required to contain the 
variably dimensioned arrays 


Routine to read 
from cards the 
first set of 
responses to be 
sorted 




Error 
return 


SRTFRP 


Figure 5.— Macro Flow Chart of Overlay (L225, 1,0) SORT 
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Figure 5.— (Concluded) 




I/O Devices of SORT 


SORT reads card input instructions (card sets 3.0 through 6.0) and reads the input 
responses from either cards (cards 4.1 through 4.4) or a magnetic file, NTAPE. 

SORT prints a summary of the items sorted and writes the items to be analyzed onto 
the magnetic file SRTFRP. 


3.3 OVERLAY (L225,2,0) - ANAL 
Purpose of ANAL 

ANAL is the second primary overlay in program L225 (TEV126). ANAL reads card 
input data describing the type of forcing function to be used in convoluting the 
frequency response functions and defining the types of output results desired. It reads 
the frequency response functions from the scratch file SRTFRP, generates the required 
time histories, prints the items requested, and optionally writes data files for plotting. 

Analytical Steps of ANAL 

ANAL performs its task in the following steps: 

1. Subroutine CARDRD is called to read card input instructions. 

2. Frequencies are read from the scratch file SRTFRP. 

3. Subroutine INITAL is called to defined execution control data. 

4. Subroutine TFRESP is called to define the transform of the forcing function. 

Steps 5 through 8 are repeated until all frequency response functions on SRTFRP have 
been processed. 

5. The real and imaginary parts of a frequency response function(s) are read from 
SRTFRP. 

6. The subroutine CONVOL is called to evaluate the convolution integral. 

7. The requested output is printed. 

8. The requested plot files are prepared. 

9. If plot files were generated, an end-of-file is written on them and then returned to 
the calling program. 

The macro flow chart for this overlay is shown in figure 6. The subroutines called are 
displayed in table 3. 
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CARDRD 


Routine to read all 
data cards 



Routine to read frequency 
array written in SORT onSRTFRP 


Routine to check the frequency 
array for monotoni cal ly 
increasing values 



Routine to define a consistent 
set of execution control 
parameters 




Figure 6.— Macro Flow Chart of Overlay (L225,2,0j ANAL 




Figure 6. -(Continued) 
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Table 3.— Routines Called by ANAL 


OVERLAY (L225, 3,0) 
PROGRAM ANAL 


CARDED 


FETAD + 
FETDEL + 

^ IRDCRD + 
MTONIC 
NAMFIL + 
READ TP + 


CONVOL 


INITAL 


r 

INVERS 

LINEAR 


( EQUAL 
LESS 
SET 


{ 


HARM 


( 


HARM 


MTONIC 
READTP + 


TFRESP 


WRTETP + 


( LINEAR 
TRANSF 


HARM 


* indicates a routine on the FORTRAN subroutine library . 

+ indicates a routine on the DYLOFLEX alternate subroutine 
library DYLIB. 


All others are local to L225 (TEV126) 



I/O Devices of ANAL 


ANAL reads card input instructions (card sets 8.0 through 11,0), frequency response 
functions from the temporary file SRTFRP, and optionally reads a forcing function 
array from the file NTFTP. 

ANAL prints the time histories ^and frequency response functions requested via card 
input. It also optionally writes the plot data files PLTRES and PLTRAN. 

o 

3.4 DATABASES 

L225 (TEV126) data bases include input and output files plus internal scratch files and 
common block storage. 

3.4.1 INPUT DATA 

Input data are from two sources: cards and magnetic files. 

Card Input Data 

For a complete description of all the card input formats, see section 6.5.1 in volume I of 
this document (User’s Guide). 


Magnetic File Input Data 

For a complete description of the magnetic file input data, see section 6.5.2 in volume I 
of this document (User’s Guide). 

3.4.2 OUTPUT DATA 

Output data may be of two types: printed and magnetic files. 

Printed Output Data 

For a complete description of the printed output data, see section 6.6.1 in volume I of 
this document (User’s Guide). 

Magnetic Files (Tape or Disk) 

For a complete description of the magnetic file output data, see section 6.6.2 in volume I 
of this document (User’s Guide). 

3.4.3 INTERNAL DATA 

Two methods are used to pass data from one portion of the program to another: a 
scratch (temporary) storage file and labeled common blocks. 
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Magnetic File (Scratch Disk File) 


The SORT program, overlay (L225,l,0), reads frequencies and frequency response 
functions from cards or tape, selects responses requested by the user, and writes them 
on the scratch file SRTFRP to be read and processed by ANAL, overlay (L225,2,0). 


SRTFRP is generated by L225 (TEV126) and released by the program when a run 
terminates normally. If a fatal error occurs, the subroutine RETURNF is not called, and 
the file will still exist after execution is terminated. This may be helpful for debugging 
purposes. 


SRTFRP is written with the subroutine WRTETP and read with the subroutine 
READTP. Both routines are described in reference 1. The contents of SRTFRP are 
displayed in table 4. 


Array Size 
NFREQ X 1 

NFREQ X 1 

NFREQ X 1 

NFREQ X 1 

NFREQ X 1 

NFREQ X 1 

NFREQ X 1 

NFREQ X 1 

NFREQ X 1 

NFREQ X 1 

NFREQ X 1 

NFREQ X 1 

NFREQ X 1 


Table 4.— Contents of SRTFRP - A Temporary Scratch File 


where ; 


|freq| 

|resr[ 

|resi|j^ 

{freq[ = 

^resrJj 

IRES1J2 

jRESRf^ = 


{RESif. = 

Inrq 

NFREQ = 

jRESRf 

jRESlfj^ 

jRESR|2 

{rESI}2 

NRQ 

End-of-File 

NRL 

1 


frequencies at which 
responses are defined 


real parts of frequency 
response i. 


imaginary parts of 
frequency response i. 


number of frequencies. 


number of generalized 
coordinate (q*s) 
frequency responses 
retained for analysis. 


number of load frequency 
responses retained for 
analysis . 


19 


I 


Common Blocks 


Table 5 displays the common blocks used in the program and the overlays in which they 
are used. 

LABELED common blocks are used for communication between the main and primary 
overlays, and for communication between routines in a primary overlay. Block names 
and contents are described in table 6. 

BLANK common is used in the primary overlay SORT as a variable-length working 
storage area. The main program of the overlay calculates the area required for arrays in 
the various subroutines and passes a first-word address and variable dimension for each 
array through the subroutine’s calling sequence. 


Table 5.— Common Blocks in Each Overlay 


COMMON 

^\BLOCK 

OVERLAY 


ERRORS 

INDUT 

MASK 

SRTFIL 

CHKPRT 

RWBUFF 


BLANK 











(L225,0,0) 
L225 VC 


X 

X 

X 

X 





(L225,1,0) 

SORT 


X 

X 


X 

X 

X 


X 

(L225, 2,0) 
ANAL 


X 

X 

X 

1 

X 


X 
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Table 6.— Contents of Common Blocks 


LABELLED COMMON NAME: ERRORS 


DESCRIPTION: Contains the number of fatal errors diagnosed 








0. 

VARIABLE 

TYPE 

DIM. 

ENG . 

NOMENCLATURE 

DESCRIPTION 




! 



NERROR 

I 

1 

! 

Number of fatal 
errors diagnosed by 

the program. 





Table 6. — (Continued) 


IKOUT 


standard input/output file names 




- — — , 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 




1 

1 

1 


The name of the card 
input file (=5) 

1 The name of the 
printed output file 
( = 6) 

The name of the punch 
file (=7) 


NOTE: 


IPFIL is not 
actually used 
inside L225. 












Table 6.— (Continued) 


LABELLED COMMON NAME: CHKPRT 


DESCRIPTION: Contains an option requesting checkout printing 


NO. VARIABLE TYPE 


ENG. 

NOMENCLATURE 



DESCRIPTION 


Option requesting 
checkout intermediate 
printing. 

In SORT 

=0 no checkout print- 
ing. " 

=1 print responses 
input to SORT. 

=2 print frequencies 
and responses in- 

-I- 4- 






LABELLED COMMON NAME:_ 
DESCRIPTION: Co ntains 

NO. VARIABLE 
1 IBl 

IB2 I I 


TYPE 


H 


RWBUFR 


R 


Table 6-— (Continued) 


RWBUFF 


a buffer array for the subroutine READTP 


DIM. 


ENG. 

NOMENCLATURE 


DESCRIPTION 


1 


Keyv;ord telling the 
routine READTP that 
the buffer size is 
word two 


in 


(=8HBUFFSIZE) 


1 

250 


Length of the follow- 
ing buffer. 

The buffer used by 
READTP to read the in^ 
put frequency and 
response matrices. 
Size of buffer limits 
the size of arrays 
read. 



Table 6.— (Concluded) 


LABELLED COMMON NAME: 


Blank Common 


DESCRIPTION : Contains t he SORT variably-dimensioned arrays 


— — — 1 

- — ^ ^ — 1 





NO, 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 

• ■ ■ ■' — 

■ 





1 

FREQ 

R 

NFREQ 


The .array of frequencies 
at which the frequency 
responses are defined. 

2 

RESR 

R 

(NFREQ, 
HRMAX ) 


The array to contain 
the real part of the 
frequency responses. 
Each column represents 
a frequency and each 
row represents a 
response . 

3 

RES I 

R 

( NFREQ , 
NRMAX) 


The array to contain 
the imaginary parts 
of the frequency 
responses. Each 
column represents a 
frequency and each row j 
represents a different 
response. 
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4.0 EXTENT OF CHECKOUT 


L225 (TEV126) was checked out with 17 different data cases using frequencies and 
frequency response functions generated by L221 (TEV156) (ref. 3). The results from 
L225 (TEV126) match those produced by an earlier version of the program.' 

One pass through SORT, the 1,0 primary overlay, sorted the response and saved them 
on a temporary scratch file SRTFRP. Then, the 2,0 primary overlay ANAL was called 
17 times, once for each data case. 

A separated deck of card input data was used to trigger most of the program’s 
diagnostics. $CHECKOUT (card 1.2) was included in the deck to allow the program to 
continue processing even though fatal errors have occurred. 

A complete list of card input data, control cards, and printed output for all test cases is 
saved on the program’s master tape and packaged with the latest program listing. 
Table 7 provides a summary of the option combinations tested. 

Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle. Washington 98124 
October 1977 


1 


See footnotes 1 and 2 in section 2.0. 
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Table 7.— Summary of Checkout Cases 









OPTIONS 

USED IN 

SAMPLE RUNS 














^ 

4 

5 

6 

7 

8 

9 

mm 

mm 

mm 

mm 

14 

mm 

mm 

17 

FUNCTion 

TYPE 

!Ii! 

■niw 

SINE WAVE 


(1- 

COS) 


THREE ONE 



NON- 

STANDARD 




IMPULSE 


Time (T) 






0.0, .82 


0.0, .82 


^91 


1^1 

1^1 

IB 

Bli 




Amplitude (AMP) 




1.0 


1.0 

1.0 




- 

- 

- 

- 

- 

- 

1.0 


Period (PERIOD) 

■HH 


HBH 

.82 

.82 

.82 

.82 

^ .82 

.82 

- 

- 

- 

- 

- , 

a . 

- 

- 

TAPE 


Source of 


- 

- 

- 

- 

- 

- 


- 


- 

- 

- 

- 

- 

YES 

- 

CAROS 


time functions 


- 

- 

- 

- 

- 

- 

- 

- 

YES 

YES 

YES 

YES 

YES 

YES 

- 

- 

0PTI0ns 


YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 


(IPTION) 

6 

6 

6 

5 

6 

6 

6 

6 

6 

1 

2 

3 

4 

7 

8 

6 



(M) 

10 

10 

10 

10 

10 

10 

10 

10 

10 

0 

6 

6 

6 

6 

6 

10 

10 


(TMAX) 

3.0 

3.0 

3.0 

0.0 

3.0 

3.0 

3.0 

3.0 

3.0 

3.0 

3.0 

3.0 

0.0 

3.0 

3.0 

3.0 

3.0 


(DELTAT) 

0.0 

0.0 

0.0 

.001 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.05 

0.05 

0.05 

0.0 

0.0 


(FHARH) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


(FMAX0) 

15.0 

15.0 

15.0 

15.0 

15.0 

15.0 

15.0 

15. G 

1 5.G 

15.0 

15,0 

15.0 

15.0 

15.0 

15.0 

15.0 

15.0 

CONJUgate 


- 

- 

- 

- 

- 

- 

- 

- 

YES 

- 

- 

- 

- 

- 

- 

- 

- 


CWPLex 

YES 

YES 

YES 

YES 

YES 

- 

- 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

INTEGrate ^ 

REAL 

- 

- 

- 

- 

- 

YES 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 


IMAGInary 

- 

- 

- 

- 

- 

- 

YES 


- 

- 


- 

- 

- 

- 

- 

- 

PRINT changes 


YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 


YES 

YES 

YES 

YES ^ 

_Y«_ 

- 

OUTPUT 

STANOard 

- 

- 


- 


1 YES 

1 YES 

- 

- 

YES 

YES 

YES 

YES 

YES 

YES " 

yes' 

- 

N0N-Standard 

YES 

YES 

YES 

j YES 

YES 

- 


YES 

YES 

- 

- 

- 

- 


- 

- 

YES 


TINE 

YES 

YES 

j 

YES 

YES 

YES 

YES 

1 YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

YES 

PRINT 


INTERpolated 

YES 

YES 

1 ~ 

- 

- 

- 

1 

1 

- 

- 

- 

- 

- 

- 

- 


YES 


ALL Transfer 

- 

- 


- 

- 

- 

1 - 

- 

YES 

- 

- 

- 

- 



- 

- 


HAXlMum 

YES 

- 

YES 

- 

YES 

- 

- 

YES 

- 


- 

- 

- 

- 



YES ' 

SELECT 




- 

YES 

- 

- 

- 

- 

YES 

- 


- 

- 


- 

■- 

- 

- 



(NTPTS) 

- 

- 

1000 

- 

- 

- 

- 

28 

- 

- 

- 

- 

- 

- 






(DELT») 

- 

- 

.003 

- 

- 

- 

- 

0.0 
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